Effects of Interband Transitions on Faraday Rotation in Metallic Nanoparticles 
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The Faraday rotation in metallic nanoparticles is considered based on a quantum model for the 
dielectric function e(u) in the presence of a DC magnetic field B. We focus on effects in e(uj) due 
to interband transitions (IBTs), which are important in the blue and ultraviolet for noble metals 
used in plasmonics. The dielectric function is found using the perturbation of the electron density 
matrix due to the optical field of incident electromagnetic radiation. The calculation is applied to 
transitions between two bands (d and p, for example) separated by a gap, as one finds in gold at 
the L-point of the Fermi surface. The result of the DC magnetic field is a shift in the effective 
optical frequency causing IBTs by i^isB/fi, where opposite signs are associated with left/right 
circular polarizations. Faraday rotation for a dilute solution of 17 nm diameter gold nanoparticles 
is measured and compared with both the IBT theory and a simpler Drude model for the bound 
electron response. Effects of the plasmon resonance mode on Faraday rotation in nanoparticles are 
also discussed. 

PACS numbers: 77.22.-d, 78.20.Ls, 78.67.-n, 78.67.Bf, 

Keywords: quantum Faraday rotation, dielectric function, interband transitions, nanoparticles, plasmons 
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I. INTRODUCTION: FARADAY ROTATION 
ENHANCEMENT AND PLASMON MODES 



There is great interest to design new materials with en- 
hanced Faraday rotationtl (FR); these mediacTEl are good 
candidates for applications like field detectors, phase 
modulators and optical isolators. Nanoparticles (NPs) 
of radius a much less than the wavelengths A of the 
electromagnetic (EM) radiation are particularly inter- 
esting possibilities,!^! because one hopes to be able to 
tune their fundamental physical properties that deter- 
mine the dielectric permittivity e(uS), which influences 
the FR signal. Further, metallic NPs or NPs withumetal- 
lic shell coatings have a surface plasmonic modeQ where 
the electron response is g*eatly enhanced, which leads 
to increased FR response!! In magnetic core NPs with 
metallic shell coatings the surface plasmonic mode in- 
teracts with a transition in the magnetic core, which is 
another process that leads to enhanced FR.Q Faraday 
rotation and cirpular dichroism in NP aggregates of var- 
ious gcomctricsEJ offers promise for control of magneto- 
optical effects. Backscattering of light in disorderedjne- 
dia may lead to enhanced Faraday rotation effects. liil 

It is simplest to use a classical Drude term to rep- 
resent approximately the dielectric response of bound 
electrons J±. or even ignore the bound electron dynamic 
response, E£l however, these approaches do not describe 
the dielectiic properties well at higher frequencies. Hui 
and Stroudtj have considered the FR response of a dilute 
suspension of small particles, with a Drude approxima- 
tion for NP dielectric function. It is the goal here to com- 
pare a classical phenomenological Drude approach with 
a quotum model appropriate for noble metals such as 
gold,oli3 where interband transitions (IBTs) take place 
from d to sp bands at the L-point £3 We consider a case 
where the NPs are sufficiently small, so that the primary 



effect of the DC magnetic field that produces Faraday 
rotation is a Zeeman splitting of the band states, rather 
than an entire series of Landau levels. EZl 

The surface plasmon frequency ui sp in a NP in the 
Rayleigh limit (a <C A) is considerably less than the bulk 
plasmon frequency lu p for the same metal. This is due 
partly to a geometrical effect, but the more significant 
reason is that interband transitions taking place above 
a gap energy greatly modify the dielectric function in 
the region of the plasmon resonance. For gold, the bulk 
plasmon, well into the ultraviolet at 138 nm, is moved 
to around 520-532 nm for the surface plasmons of NPsJ3 
But to get the correct description, the interband tran- 
sitions must be taken into account for describing e(uj). 
Here, we include IBTs for the bound electron contribu- 
tion to e(w), in the presence of a DC magnetic field, so 
that the Faraday rotation properties can be described 

Faraday rotation is a magneto-optical phenomenon!! 
that measures the fundamental electronic, optical and 
magnetic response-of a dielectric medium. It is similar 
to optical rotationEJ except that FR requires an applied 
magnetic field. The Faraday rotation is the change in 
the polarization of an EM wave as it propagates through 
some medium parallel to the axis of a quasi-static mag- 
netic field. The interaction between the DC field and the 
charges in the medium leads to different speeds of propa- 
gation and different wave vectors kn and k^ for the right 
and left circular polarization components of EM waves, 
leading to the net rotation of initially linearly polarized 
waves. 

The basic parameter to describe the degree of FR is 
the Verdet factor, v, which is the rotation angle of the 
polarization per unit propagation length z per applied 
magnetic field B: 
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We consider a composite medium of gold NPs in water. 
At small volume fraction of gold, f s -C 1, the Verdet 
factor is linear in the volume fraction. Generally, the 
rotation ip is also linear in magnetic field at lew enough 
fields (although materials of other symmetry^ can ex- 
hibit quadratic dependence on B). Then v does not de- 
pend on B, and we consider only this regime. For dilute 
composites, the Verdet factor per volume fraction is a 
better quantity for consideration, defined as 

fs Bzf s V ' 

The Faraday rotation results from the phase difference 
of the polarizations, 

<p = ^Re{k R - k L }z , (3) 

where the respective wave vectors for propagation of the 
two circular polarizations are determined by relative di- 
electric functions en and ex, 

l = —JWr^E ■ (4) 
c v 

with ji being the relative magnetic permeability of the 
medium and c being the speed of light in vacuum. The 
theoretical description of the FR signal in some metallic 
or metallic-shell NPs, as a result, requires an accurate 
description of the dielectric function e(w) for the metal, 
including the presence of the DC magnetic field. Espe- 
cially, it is important to have a reliable description of 
e(w) at shorter wavelengths, into the ultraviolet, below 
the plasmon frequency for the NPs. Indeed, in order for 
the theory to correctly predict the plasmon frequency 
Lo sp , requires knowing e(w) all the way into the ultravio- 
let. Thus, the goal here is to get an accurate theory for 
e(uj) with the DC magnetic field present, that includes 
interband transitions for the bound electrons as well as 
the usual plasmon response of the free electron gas. 

Below we begin by describing the NP synthesis and 
a description of the experimental measurement of the 
Faraday rotation. Then we will continue by summarizing 
the basic relations among the dielectric functions e(o>), 
ex, (a;), and ex^(w) and the FR response. (All e are under- 
stood to be relative dielectric functions, leading to index 
of refraction, n = i/fle.) Next, a simple classical model 
for e(u>) using bound electrons and based on the Drude 
model is described briefly, for comparison with the quan- 
tum calculation for the IBTs. In that model the effect 
of IBTs is approximated by the response function of a 
set of bound electrons, with some binding frequency ljq. 
Part of the motivation for the quantum calculation of the 
IBTs is to determine the validity of this simpler classical 
model. 

For the quantum effects of IPTs, we adopt the axt 
proach used by Boswarva et a/.E3 and also by Adlcio 
of finding the perturbations of the electron density ma- 
trix that are caused by only the electric field of the EM 



waves. The optical magnetic field is ignored. However, 
the DC magnetic field enters because it shifts the band 
states. This is a simple Zeeman shift; for nanometer- 
sized systems there is no sense to Landau levels that 
were used by Boswarva_e< ai, and also in theory devel- 
oped by Halpern et al.jEn due to the geometric confine- 
ment (the NP radii are smaller than the Landau radius 

r ° = \J iff' f° r a PPh e d magnetic induction B). Follow- 
ing through the calculation, there is a contribution to 
e(w) that requires summing over IBTs with a range of 
energies. Those integrals are evaluated in two different 
models: A three-dimensional (3D) band model and a ID 
band model. The results are presented both with the 
presence of a phenomenological electron damping con- 
stant 7 and in a limit that this damping goes to zero. 
The results are also compared with earlier calculations 
of the IBT contribution by Inouye et alS3 and by Scaf- 
fardi and Tochojlj that did not include the DC magnetic 
field. 

The net e(u>) includes contributions from both the 
bound (IBTs) and free (electron gas) electrons. We ap- 
ply the results to calculate the scaled Verdet factor T for 
a dilute solution of solid metallic NPs, using the parame- 
ters for gold in water. The effects of the dilution are pap, 
sidered most simply by the Maxwell-Garnett theory,E3lHJ 
assuming that the NPs do not aggregate. We find that 
even with the quantum IBTs included for bound elec- 
trons, the experimentally measured FR signal is about 
10 x stronger than that predicted by the theory. 

We will conclude with comments on the applicability 
of the results for other systems with plasmonic enhance- 
ments of dielectric responses. 

II. SYNTHESIS OF GOLD NPS AND FARADAY 
ROTATION MEASUREMENTS 

The large gold nanoparticles are prepared from the re- 
duction of HAuCL; solution by-sodium citrate solution as 
described by Turkevich et alSB Briefly, 5 mg of HAuCL 
and 50 mg of sodium citrate are dissolved in 95 ml and 5 
ml of doubly distilled water, respectively. The HAuCLi 
solution is heated to about 70° C and the sodium citrate 
solution is added, vigorously stirring the solution for 50 
minutes. The color of the solution gradually changes 
from faint pink to wine red. The resulting large gold 
nanoparticles have size 17 ± 3 nm. Assuming 100% re- 
duction of the gold into NPs, the upper limit of volume 
fraction of gold in the solution is f s = 1.50 x 10~ 6 . By 
analyzing the extinction coefficient of the solution by the 
techniques in Ref. the actual volume fraction is esti- 
mated to be f s = 1.23 x 10~ 6 . 

The Faraday rotation spectrum of NPs in water solu- 
tion was measured with the help of a home built pulsed 
magnet. The magnet consists of a helical coil machined 
from a copper beryllium block and electroplated with 
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silver. The pulsed current to the coil is provided via 
a simple RLC circuit. The capacitor bank of 77.3 /iF 
from Maxwell Laboratories is charged by a power sup- 
ply/charger of Lumina Power, Inc. The power supply 
uses 100-240 V AC-50/60 Hz input and output of 10 kV 
at 500 J/s in continuous operation. The charge from the 
capacitor bank is discharged into the coil via a high volt- 
age trigger spark gap. The current is monitored in the 
circuit via a Rogowski coil, which measures the current 
derivative. 

The Faraday rotation of the nanoparticle solutions is 
measured in a plastic cell placed in the coil. A flash 
light source is triggered along with the pulsed magnet 
that allows the synchronization of the magnet with the 
optical measurement. The duration of the light pulse 
is 1.4 /is, while the duration of the magnetic pulse is 
~ 50 /is, which allows that during the optical measure- 
ment the magnetic field is relatively constant. In front 
of the flash light source a polarizer is placed to produce 
polarized light for the Faraday measurement. The po- 
larized light passes through the sample containing the 
nanoparticle solutions. The light leaving the optical cell 
passes through another polarizer that is set to 45 degree 
with respect the first polarizer. The light then enters a 
fiber optic spectrometer, which is also synchronized with 
the pulsed magnet and the light source. The Faraday ro- 
tation is calculated from the intensity change in the spec- 
trum before and after the magnetic pulse. The magnetic 
field and the Faraday rotation setup are calibrated with 
water placed into the optical cell. The measurements are 
taken at 4.2 tesla magnetic fields. 



It is useful to realize another way to get to e, via aver- 
aging of the microscopic currents, i.e., those caused by 
the optical fields. The dielectric medium under study 
has current density J, which combines with the vacuum 
displacement current. In this view the Ampere/Maxwell 
Law is 



<9E 

V xH = J + e - 

ot 



(8) 



All the effects of the medium are contained in J. This 
must be equivalent to the alternative viewpoint that the 
currents are represented instead by a dielectric function, 



V x H = 



<9D 

~dt 



(9) 



Considered at the frequency of the EM radiation with 
time derivatives d/dt — > — iu>, these alternate views give 



J = — icuEQ (e — 1) • E = — iuix ' e oE 



(10) 



Thus, an averaging of the microscopic currents will also 
lead to the susceptibility and dielectric tensors. 

We assume that the DC magnetic field B is along the 
z-direction, the same as the propagation direction of the 
EM waves, with wave vector k = kz. Then the elec- 
tric field in the waves has only xy components; only the 
transverse part of the dielectric tensor isjieeded. In this 
situation it has the following symmetryEZI 



£-xx ^xy 
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(11) 



III. THEORY: DIELECTRIC POLARIZATION, 
CURRENTS, e{oj) AND FARADAY ROTATION 

We consider EM radiation at frequency ui with the 
electric field E(i) ~ e~ ZUJt , incident on a material particle 
(an individual NP) much smaller then the wavelength 
(the Rayleigh limit). Then the field E is taken as uniform 
inside the sample. The dielectric properties are based on 
the averaged dipole moment of the electrons of charge e, 
d = er. For n — N/V electrons per unit volume, the 
electric polarization can be expressed as 



X-e E 



(5) 



where eo is the permittivity of vacuum and x is t ne sus- 
ceptibility tensor that is to be found. The dielectric func- 
tion considered as a tensor e is defined via the electric 
displacement D = eqc • E or 



D = e E + P , 

from which the usual definition results, 

e= 1+X ■ 



(6) 



(7) 



The off-diagonal elements are determined by the DC 
magnetic field; they vanish when B = 0. The vari- 



able £ x 



-ie a 



is convenient later; it is real in the 



absence of electron damping. The EM waves that prop- 
agate without any change in polarization are those with 
polarization vectors that are eigenvectors of e. These 
eigenstates are the usual states of right and left circu- 
lar polarization. Thus, solving the eigenvector problem, 
e • iii = ejUj, i = 1,2, with eigenvalues q and eigen- 
vectors Ui, one finds the right circular polarization state 
(negative helicity) with E y = —iE x : 



' &XI1 



1 

V2 V 



p-iy) , (12) 
and the left circular polarization state (positive helicity) 



with E y = +iE x 



(13) 



Each mode has a different wave vector for propagation, 
according to expression (|). Then starting fr om a lin- 
early polarized wave at position z = 0, its right and 
left circular components get out of phase by the time it 
travels to position z, leading to the rotation of the polar- 
ization through the angle <p given in expression (0). One 
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might also mention, that in general, the dielectric tensor 
elements are complex, then there is also a change in el- 
lipticity X of the polarization, given from the imaginary 
part, 



X = ilm (kfi — fci,) z 



(14) 



The two effects of Faraday rotation and change in ellip- 
ticity (tan X = ratio of minor to major axis of the ellipse 
swept out by the electric, vector) can be combined into 
one complex parameter J23 



ip = <p + iX = | (k R — k L ) z 



(15) 



Usually these effects are extremely small and close to 
linear in B. Then there is only a tiny difference in k R 
and fcj,, which gives to a very good approximation, the 
complex relation, 



ip = ip + iX « 2~\l- — £xy 



(16) 



This emphasizes how the components of e are needed to 
describe the changes in the optical polarization. 

From the experimental perspective, the measurement 
of the absorption (or, attenuation) coefficient a is at least 
one technique that sets a relative scale for the FR. It is 
given from 



a = 2Im{fc c ff} = 2— 1m{W{ie e s} 
c 



(17) 



This could use either cr or or their average for the 
effective dielectric function e e ff of the medium, as this 
expression does not involve their difference, which is ex- 
tremely small. Thus, measurements of a serve to set 
some unknown fitting parameters, when needed. 



IV. CLASSICAL PHENOMENOLOGICAL 
MODEL FOR e(u) (DRUDE MODEL) 

In this section the electron motion is assumed to be 
classical. An electron of bare mass m and charge e = 
-1.602 x 1CT 19 C has some trajectory r(t) = (x(t),y(t)) 
in response to all forces acting on it, and the averaging 
of its induced electric dipole moment d = er lead to the 
dielectric function. 

To include the effect of the constant B on e it is as- 
sumed that there are two primary contributions to the 
dielectric response. The first is the contribution of free 
electrons with number density n, and some damping pa- 
rameter 7 P , that leads to the usual plasmon response 
with a plasma frequency lo 2 = ne 2 /me - The second is a 
contribution due to bound electrons, with some binding 
frequency loq and another damping parameter 70. The 
contribution o£,bound electrons is essential to describe 
e(w) correctlyta in NPs. 



Any electron, whether free or bound, is acted on as 
well by the electric force from the optical held, and the 
Lorentz force from the DC magnetic held. The force due 
to the optical magnetic held can be ignored in lowest or- 
der. In this Drude approximation the equation of motion 
of a bound electron iscJ 



m n r = eE + er x B 



m„w^r - 



m o 7 r 



(18) 



Under the assumption of e %ut time dependence of the 
optical field E, which is the source field, this is 

[to g (wq — lo 2 — iw7o) — iweBx] r = eE . (19) 

In terms of the components this is a matrix relation, 



(20) 



LO 2 — LO 2 — lW7o +ILOLOB 




X 


e 


E x 


—ilolob lo 2 ~ lo 2 — iujjo 




. y _ 


m a 


_ E y 



where the cyclotron frequency with B along z is 



LO B 



eB 
m„ 



(21) 



The matrix f2 2 on the LHS of ( po| ) has the same kind 
of symmetry as that of e in (jll]), because the diagonal 
elements are equal and the imaginary off-diagonal ele- 
ments differ only in sign. This means fl 2 has the same 
eigenvectors, which are the right and left circular polar- 
ization states. Based on its structure, the eigenvalues 
Slfj and Vl 2 L of f2 2 are easy to read out. For right circular 
polarization, 



n 2 



lo" — iuj'-fo + lolob 



(22) 



For left circular polarization, the last term (from the 
off-diagonal element) has the opposite sign, 



ill 



lLO"fo — tOLOB 



(23) 



The effect of the DC magnetic field appears only in the 
last factor. These two eigenvalues can be combined into 
a single convenient expression in terms of the helicity 
v = — 1 for right circular polarization and v = +1 for 
left circular polarization: 



n. 



ILO^fo — VLOLOB 



(24) 



The helicity is the projection of the photon intrinsic an- 
gular momentum L on the direction of propagation (k 
or 2). In this expression it multiplies the magnetic held 
component along the same axis. Any physical differences 
for right and left circular polarizations will become inter- 
changed if the direction of the magnetic field is reversed. 
In terms of a vector v = L/H, the last factor in these 
eigenvalues could be written most generally as lov ■ lob- 
If the electric field contains only one of the circular 
polarizations, i.e., E = E v u u , the response r also will be 
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proportional to the same eigenvector. Then the solution 
for the electron position is very simple, 

eE v 



The factor e/mH, 2 gives the size of the response for this 
polarization. Here we see the fundamental physical dif- 
ference between the polarizations. One of the polariza- 
tions causes a larger circular motion of the electrons than 
the other polarization. Which one actually is larger de- 
pends on the relation between luq and uj. This difference 
leads to a corresponding difference in the dielectric ef- 
fects. 

Based on this position response, it is then easy to find 
the effective dielectric functions for the two circular po- 
larizations, using d = er. The susceptibility due to these 
electrons is \ = n (d)/eo-E. The result can be summa- 
rized in a single formula, 



(25) 



(26) 



where v = —1/ + 1 refers to R/L polarizations, respec- 
tively. This applies to separately, the contribution from 
the bound electrons, or, the contribution of the free elec- 
trons, using appropriate parameters in each case. 

Look at this another way. An arbitrary electric field 
can be expressed either as E = E x x + E y y or as E = 
ErUr + ElUl, where 



Er 



^(E x + iE y ), E X = ±(E R + E L ), (27) 



^ (E x - iE y ) , E y = ^ {E R - E L ) . (28) 

One can combine the right and left solutions and get the 
general solution for any electric field, in diagonal form: 



eE 



R 



o2" U « 

•m \rr 



eE L . 

7^2 U L 

m \r. 



(29) 

"R iiwl 

Alternatively, this can be written in Cartesian compo- 
nents, 



x 



1 e 

y/2m 
1 ie 



Er El 



n 2 

iL R 
Er 

Vr 



El 



Simplifying, or inverting the matrix equation (|20| 
to the general electron motion, 



X 


_ e/m 


. y _ 


^R^L . 



-IUJUJb 



-IUJUJB 

uj 2 — iujjo 



(30) 

(31) 
leads 

E x 

En 



(32) 

Multiplied by e, the relation shows the polarizability ma- 
trix of the electron. This expression leads to the suscep- 
tibility tensor, 

2 / r 9 9 

ne /m eo oJq — oj — tujjo —iujlub 

-iujuib Wq — uj 2 + icu^fo 



X 



^L^R 



One can see this is consistent with (|26|), because its eigen- 
values are \R = Xxx - iXxy and xl = Xxx + iXxy, which 
agrees exactly with (E6h. 



A. Combination of free and bound electron 
responses 

Now to use this to describe a metal such as gold, we 
assume first there is some density of free electrons n, 
with a bulk plasma frequency uj 2 = ne 2 /m eo, a damp- 
ing 7 P and a zero binding frequency. In addition, there 
is some other density no of bound electrons, leading to 
an effective weight g 2 , = n e 2 /m eo, with an associated 
damping 70 and binding frequency ujq. The net dielec- 
tric function is the sum of the two contributions to x- 
In terms of the polarization states v = ±1, the dielectric 
function is taken as 

/ \ 1 u l 9o 

e v {w) = 1 7. ; 5 5 ; 

U!" + lUlJp + VLUUJB OJ — LUq + tOJjQ + VUJUJB 

,( 34 ) 

The first two terms are the usual ones for describing a 
free electron gas. The last term uses the single resonance 
to approximate the effects of bound electrons. Both in- 
clude the DC magnetic field implicitly in the cyclotron 
frequency, ojb- The ease with which the magnetic field 
is included in the bound electron response is the main 
advantage of this model. 

One can produce the Cartesian elements of e, for in- 
stance, using ([[2]) and (|l3|), by the combinations of these 
eigenvalues: 



^X X 

e xy 



y — 2 ( e R + 6 l) , 

E yx = i£ xy = § (e_R — £l) 



(35) 
(36) 



(33) 



B. Maxwell-Garnett averaging for dilute solutions 

The medium of interest here is actually a dilute so- 
lution of NPs at a volume fraction / s < 1 in a host 
liquid, which we take as water, with its host dielectric 
constant eh — 1.777. The NPs are considered the scat- 
terers with dielectric function e s For comparison with 
experiment, the effective dielectric function e c ff of the 
solution is required. The theory for calculating the ef- 
• fective dielectric function depends somewhat on the as- 
sumption of how the particles are dispersed in the liq- 
uid. In the simplest approximation, they are assumed 
to be randomly dispersed and not sticking .to each other. 
In the Maxwell-Garnett (MG) theoryoo one finds the 
volume-averaged electric field and the volume-averaged 
polarization response to that field, from which e e a is de- 
termined. The MG theory is known to-apply well even 
in the presence of multiple-scattering£9 These volume 



averages are 



E av = / S E S + (1 — /s)E/j , 
P av = f s P s + (1 - /,)P h , 



(37) 
(38) 



where s and h refer to the values in the scatterer and 
the host, respectively. For spherical scatterers exposed 
to asymptotic field E/, in the host, the Clausius-Mosotti 
equation gives the internal fields, 



E s 



3e/, 



2e h 



Then with polarization 
one finds the average 

P 

Eeff = 1 



P a = (e a -l)e E s . (39) 
(eh — l)eoEft in the host, 

1 + 2/?/ 



£oE a 



eh 



l-(3 f 



(40) 



which involves the scaled volume fraction (/ s is the frac- 
tion of volume occupied by NPs in the solution), 



Pf = f S 



2e h 



(41) 



This MG averaging procedure for composite systems is 
usually summarized by the equivalent relation, 



e e ff + 2eh 



- = fs 



e s - e_h 
e s + 2th 



(42) 



Expression (40) can be applied separately to the left and 
right circular polarization states, then leading to an ef- 
fective dielectric function for each, that will then give 
the Faraday rotation (0) for a dilute solution. 



C. Classical model parameters for gold 
nanoparticles 

Based on the work in Ref. ^, the parameters needed 
for this classical model were found by fitting it to the ab- 
sorption measured experimentally with B = 0, for a di- 
lute solution of 17 nm diameter gold NPs in water. That 
fitting is based on using the effective dielectric function 
e e ff from the MG theory, to give the absorption in the 
solution, according to expression (|T7|). 

For this classical Drude model, based on the elec- 
tron number density, and using effective mass equal to 
the bare electron mass, the bulk plasma frequency is 
ojp = 1.37 x 10 16 rad/s, which corresponds to \ p = 
2irc/ujp — 138.5 nm. The damping of the free electrons in 
NPs can have an intrinsic term and a surface scattering 
term. Thus a size-dependent damping factor is included, 
according to the combination of these processes,!^ 



7 P 



Vf 

T 



(43) 
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17 nm gold particles in water, 
absorption over 1 cm distance 

free electrons: 

CD = 1.405 x 10 16 rad/s (X - 134.1 nm) 
Y p = 1.118 x to' 4 rad/s (T. = 8.94 fs) 

bound electrons, Drude model: 
g = 3.698 x 10 15 rad/s 
ffl = 3.734 x lO 15 rad/s (504.4 nm) 
y = 6.084 x 10 14 rad/s (1.64 fs) 



experimental data 

fit: bound + free electrons 
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FIG. 1: (Color online) Fitting of the absorption of 17 nm 
diameter gold particles in water solution, according to the 
Drude model for the bound electrons. Parameters indicated 
are used to get a good fit to the absorption peak near 522 
nm. The fitted volume fraction of gold is f s = 3.36 x 10 -6 . 
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FIG. 2: (Color online) The bound electron contribution to 
the permittivity, from the second term of Eq. (|3l|) , according 
to the Drude model for the bound electrons, using the param- 
eters of Figure |lj The real part of e(ai) becomes negative for 
frequencies above ojo (wavelength 504 nm), which is a defect 
of this model. 



where r « 9.1 fs is the intrinsic scattering time, vf = 
1.40 x 10 6 m/s is the Fermi velocity, and d is the thickness 
of the gold. This thickness could be the diameter for 
solid spherical particles, or, the thickness of a shell for 
core/shell particles. We discuss data for gold particles of 
average diameter 17 nm; the prediction for their effective 
damping is then 7 P = 1.92 x 10 14 rad/s, or, a time scale 
t p = 7- 1 = 5.20 fs. 

The Drude theory was fitted to experimental data for 
absorption through a 1 cm path of water solution of 
gold particles with average diameter of 17 nm. The 
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FIG. 3: (Color online) Faraday rotation for 17 nm diame- 
ter gold NPs, from experiment, and according to the Drude 
model for the bound electrons, using the Drude fitting pa- 
rameters of Figure jjj. The experimental results have been 
scaled by 1/100 to allow them to be plotted together with 
the Drude theory. The theory result was obtained with the 
MG effective medium approach. Faraday rotation angle ip 
and ellipticity angle X have been scaled by B,z,f s , to give 

■$> = r + iz. 



fitting parameters were chosen to get a good descrip- 
tion of the absorption peak present near 522 nm, at- 
tributed to surface plasmon response. A good descrip- 
tion can be obtained while also allowing the volume frac- 
tion and free electron parameters uj p and j p to vary, 
see Figure [I]. The contribution from the bound elec- 
trons can be represented approximately using the am- 
plitude parameter g = 3.70 x 10 15 rad/s, the binding 
frequency ujq = 3.73 x 10 15 rad/s (wavelength 504 nm), 
and damping frequency 70 = 6.08 x 10 15 rad/s, which 
corresponds to a damping time of r = I/70 ~ 1-64 fs. 
To get this good fit to the peak, the free electrons are 
at the same time represented using plasma frequency 
uj p = 1.40 x 10 16 rad/s, equivalent to X p = 134.1 nm, 
and a damping -f p = 1.118 x 10 14 rad/s, correspond- 
ing to the damping time t p = 8.94 fs. These are slightly 
different than the accepted bulk values, however, we con- 
sider them here only as a model that fits accurately the 
absorption peak. 

From these fitted dielectric parameters, the theoretical 
Faraday rotation response can be obtained. Results for 
the Faraday rotation and ellipticity for 17 nm gold NPs 
in solution are shown in Fig . [}| The complex rotation 
angle ip is found from Eq. [Lfj, together with applying the 
Maxwell-Garnett procedure for the composite medium, 
Eq. [lo], for the effective dielectric function of the com- 
posite solution. We have scaled the rotation angle tp and 
ellipticity X by the product of path length z, magnetic 
field B and gold volume fraction f s , to remove the lin- 
ear dependence on these quantities. Thus we define the 



complex rotation angle scaled by volume fraction, 
* = T + iZ = 1>/f, . 



(44) 



Then, T = tp/(Bzf s ) is the Verdet factor per unit vol- 
ume fraction, and Z = X/(Bzf s ) is a corresponding 
ellipticity factor per unit volume fraction. Then the re- 
sults for T and Z do not depend on B, z, or f s in the 
linear regime. The experimental data for T are also dis- 
played in Fig. ||, scaled down by a factor of 1/100 in 
order to be shown together with the theory. 

One sees that the model predicts a negative peak in 
the Faraday rotation near 525 nm, apparently associ- 
ated with the plasmon resonance (see ^8] for the dis- 
tinction between positive and negative rotation angles). 
The experimental data have a similar negative peak in 
the same region, although its magnitude is significantly 
larger than this theory predicts. The theory has a wider 
positive peak around 580 nm and a long tail at longer 
wavelengths, but this positive peak is rather weak in 
the experimental data. For the ellipticity, the main fea- 
ture predicted is a positive peak around 540 nm, slightly 
above the plasmon wavelength, together with its associ- 
ated negative peak and long tail at shorter wavelengths. 
Unfortunately, the model exhibits an artifact at shorter 
wavelengths: both tp and X tend to increase greatly at 
short wavelengths in an unphysical behavior. This is due 
to the fact that classical Drude model cannot correctly 
describe the bound electron response at higher frequen- 
cies. 

This model is an approximate way to include the ef- 
fect of B on classical bound electrons, however, it should 
be replaced by the more complete calculation using the 
quantum interband transitions presented later. It gives 
a reasonable fit to the absorption curve from 900 nm 
down to 400 nm, however, below that wavelength it pre- 
dicts much more absorption than actually takes place. 
Also, this Drude description of the bound electrons can- 
not accurately describe the response in the wavelengths 
350 - 500 nm. This model does not require any back- 
ground (i.e., high-frequency) dielectric function f M ~ 10, 
as has been applied-in other studies to mimic the effect 
of bound electrons. U Even so, the fit to the absorption 
peak due to the SP mode is very good, while the corre- 
sponding negative FR peak due to the SP mode of 17 
nm gold particles is about 100 times stronger than the 
theory predicts. 



V. QUANTUM DESCRIPTION OF e{u) VIA 
PERTURBATION OF THE DENSITY MATRIX 

In this Section we consider the quantum calculation 
of the effects due to bound electrons, which is taken into 
account by finding contributions to e(w) due to interband 
transitions, in the presence of the DC magnetic field. 
The electrons are considered non-interacting. 
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The single-electron Hamiltonian is taken as 



H 



1 



2m D 



p-eA(f,i) +e(j)(r,t) + U(r) , (45) 



where the charge is e, <j> and A are the scalar and vec- 
tor potentials of the EM fields, and U is the periodic 
potential of the lattice. The canonical momentum oper- 
ator for the electron is p = — iftV. The more physical 
momentum is the kinetic momentum operator, 



eA 



(46) 



because it is the square of this operator that deter- 
mines the energy. We take the scalar potential as zero 
(Coulomb gauge). The vector potential A includes a 
term for the DC magnetic field, Ao and a term for the 
AC optical field Ai. The optical field is treated as a 
classical non-quantized field, that oscillates as e" lwt . 

The electron bands, unperturbed by optical fields, 
come from the solution of a Hamiltonian with the ki- 
netic energy, the lattice periodic potential, and the DC 
magnetic field, 



^ = i(p- eAo ) 2 + ^' 



(47) 



At weak enough DC magnetic field, the quadratic term 
in Ao can be dropped, and the effect of the cross term 
with p is the orbital Zeeman splitting, 



H 



2m, 



(48) 



The magnetic dipole moment due to the orbital angular 
momentum is 



2m 



(49) 



As the electron charge is negative, fl points opposite to 
L. The component of jl along the magnetic field is m/is, 
where [is ~ eh/2m is the (negative) Bohr magneton, 
and m = mi is the magnetic quantum number. The 
Zeeman splitting is an energy shift AE — —mpisB = 
— ^mhwB, where both \xb and ojb can be negative, due 
to the negative electron charge. The states of this Hamil- 
tonian are some electron band states, including any Zee- 
man shifts, 



H \klm) = Ew m \\dm) . 



(50) 



The band states, labeled by wave vector k and angular 
indexes I, m, have wave functions 

tpkim(r) = (r\Mm) = -^=e <k ' r M Wm (r) . (51) 

These can be considered the original states of an unper- 
turbed problem. The optical field is the perturbation 



on these states, whose effect is studied using the density 
matrix approach. 

Because we consider states in NPs, the DC magnetic 
field only produces Zeeman shifts, rather than Landau 
levels. Due to the geometrical confinement, there is no 
sense to Landau levels that would have extended wave 
functions much larger than the size of the particles. For 
instance, at a small applied field strength B = 0.1 tesla, 
the length scale of the Landau levels is the Landau ra- 
dius, tq = ^/|§ = 115 nm. This is much larger than the 

radius of the nanoparticles under consideration, typically 
from 5-10 nm. The Landau wave functions do not fit 
into the NPs at this field strength, giving a non-bulk sit- 
uation. The degeneracy of Landau levels is on the order 
of (R/ro) 2 , where R is the system radius. At B = 0.1 
tesla, the degeneracy is about (8.5/115) 2 w 0.0055, how- 
ever, this fractional value is not meaningful. For the 
larger magnetic field B = 4.2 T, the Landau radius is 
reduced to r = 17.7 nm. This is still somewhat larger 
that the NP radius of 8.5 nm, and the degeneracy is 
about (8. 5/17. 7) 2 w 0.23, still significantly less than 1, 
so the theory should be applicable. 

These considerations show that the Landau levels are 
the incorrect solutions in a confined geometry. When one 
looks more carefully at how to arrive at the quantum so- 
lution, the radial wave functions should go to zero at the 
boundary of the NP (for bound electrons). For spheri- 
cal particles, that radial dependence would be described 
by spherical Bessel functions, ji(kr), with discrete al- 
lowed k, and angular dependence described by spherical 
harmonics for a chosen angular momentum, l,m. We 
consider a quasi-bulk approximation, where the discrete 
k are assumed to be close enough together to be reason- 
ably described by electron bands. 



A. The density operator p 

Statistically, the band states are populated according 
to a Fermi-Dirac distribution for the given temperature, 
when the system is in equilibrium. The density operator 
p is a way to introduce this population into the QM 
problem and provide for mixed states. Of course, once 
the optical field is turned on, a new equilibrium can be 
established and the density operator can change. Its 
basic definition for an equilibrium situation, in terms of 
the state probabilities Wi is 



(52) 



For the equilibrium distribution, the weights are taken 
as proportional to the Fermi-Dirac occupation numbers, 



w i = jjf (E i ), f (E i ) = 



1 



e P(Et-E F ) + l ' 



(53) 
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where /3 is the inverse temperature and Ep is the Fermi 
level. This is a density matrix normalized to one for iV 
electrons in the system. The time derivative of a general 
p follows from the quantum Liouville equation: 



at = m [H ' p] 



(54) 



To apply this, we consider the leading perturbation 
term in the Hamiltonian, which is the electric force from 
the optical field (Ai). The optical magnetic force is ig- 
nored. Then the perturbation is described by the Hamil- 
tonian 



Hi = - 



rrin V 



eA ( 



(55) 



Now the total density operator is assumed to be a sum 
of the equilibrium operator plus some change caused by 
the perturbation: 



P = Po + Pi ■ 



(56) 



As the total Hamiltonian also is a sum of unper- 
turbed and perturbation parts, we can use the fact 
that [Ho, pa] — 0, and ignore the small nonlinear term 
[Hi, pi] ~ 0, then the equation of motion for the pertur- 
bation is 



lh ~gf ~ [ H 0,Pl\ 



[Huh] ■ 



(57) 



Now assume expansions of Hi and of p\ in the unper- 
turbed basis states, 



Hi=^2\f){f\Hi\i){i\, pi = Y,c fi \m 



(58) 



if 



if 



The constants c/j = (f\pi\i) are just the matrix elements 
of pi, in the Hq basis states. After evaluation of the 
commutators, and assuming e~ lu)t time dependence for 
pi, the constants cu are found, and the change in the 
density operator is found to be 



{ Wi - ws )\f){f\Hi\i)d\ 



(59) 



A small imaginary part 7 has been added to the fre- 
quency to effect the turning on of the perturbation. This 
constant can be considered a phenomenological damping 
constant, or, it can be let to go to zero if the results with- 
out damping are of interest. This [expression has been 
used in various, problems by AdlerEi and in the thesis 
of M. Prange£3 Ei and Ef are energies of two states of 
the unperturbed Hamiltonian. One can think that the 
expression involves transitions between pairs of states. 
Obviously the oscillatory time behavior of the perturba- 
tion Hamiltonian must be reflected in a similar behavior 
in this part of the density matrix. Thus, we are inter- 
ested only in the response in the density matrix at the 
same frequency as the perturbation. 



B. Thermal and volume averages 

To find the dielectric function, statistical averages of 
the polarization or the current density are necessary. 
This can be done by first defining a local quantum op- 
erator, p e (r) for the one-electron charge density, 

p e (r)=e|r)(r| , (60) 

and another, j(r), for the one-electron current density, 

i(r) = |{|r)(r|v + v|r)(r|} . (61) 

The current density operator is defined in terms of the 
electron velocity, 



v = — = — ( p - eA ) . 

m m 



(62) 



The statistically averaged values of these operators are 
found from the trace with the density operator, 



Pe (r) = (p e (r)) = Tv{p Pe (r)} 
j(r) = (i(r))=Tr{pi(r)} ■ 



(63) 
(64) 



In a pure state \ip), with density operator p = \ip}(ip\, 
these produce the usual expressions for the quantum 
charge and current densities at point r, 



p e {x) = e|V(r)| 2 , 
j(r) = Re {-0*(r)evV>(r)} 



(65) 
(66) 



We can also define the local polarization operator using 
the electron position, 

d(r) = er|r)(r| , (67) 

which is statistically averaged by the same procedure, 

d(r) = Tr|pd(r)| . (68) 

At some point in the calculation the volume averages 
are desired, to describe e(uj) for the whole sample. These 
are obtained from the usual definition, say, for the charge 
density (the overbar indicates volume average) due to N 
electrons in a volume V, 

Pe = yJ d 3 rp e (r) = ^J2 W *J d 3 r\Mr)\ 2 ■ (69) 

i 

The individual electron states ipi are unit normalized and 
the probabilities sum to one. This recovers an obvious 
result, 



N 

P = N P e = yTr{pe} = ne 



(70) 
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For the current density and the electric polarization, 
there are similar expressions, 



J 

P = 



Nd 



d 3 r j(r) = Tr {pnev} 



N 
V 

— / cfVd(r) = Tr{pner} 



(71) 
(72) 



For the most part, we use the averaging of the operator, 
V = ner , (73) 
to determine the volume-averaged polarization response. 

C. Averaged electric polarization response 

The perturbation oscillates at frequency u and there- 
fore we get the response in the electric polarization P 
that oscillates at the same frequency, using only the 
change pi in the density matrix, 



P = Tr 



{piV]=neJ2(f\pi\i)mf) 

if 



(74) 



This requires matrix elements of the position operator 
in the unperturbed basis states. Those can be obtained 
from the equation of motion in the unperturbed system, 



ihr = ihv = [r, Hq] 



(75) 



Then the needed matrix elements can be expressed using 
the velocity, 



(i\r\f) 



ih 



(Ef - Ek) 



(*|v|/) 



(76) 



J 



The optical electric field is E = —(dA±/dt) = i(oj + 
i"f) Ai, including the turning on of the perturbation. The 
perturbation can be expressed now as 



Hi = -eAi • v = 



i(u> + ry) 



E • v , 



(77) 



The operator v need include only the DC vector poten- 
tial, Ao. Then the matrix elements of both p and Hi 
come from the velocity. The result for the averaged elec- 
tric polarization is expressed as 



(w + 17) 2-i [ft(w + %i) + E t - E f ] (Ei - E f ) ' 

(78) 

As E is assumed to oscillate at frequency ui, this is indeed 
the response oscillating at that same frequency. The 
damping 7 is necessary so that an appropriate limit gives 
the classical damped responses found earlier. 

With the transition frequencies given by 



Huiif — Ei — Ef , 
the susceptibility components that result from (p 



Xab — 



e h(uj + 17) 



E 

if 



(Wj - Wf)(i\v a \f)(f\v b \i) 

LUif (w + ij + Uif) 



(79) 
are 



(80) 



We apply this to find only the effects from interband 
transitions. The free electron response in ( |34|) is still 
applied for the quantum model. 



The result can be symmetrized by labeling some states as occupied states (o) and the rest as unoccupied (u) . All 
terms correspond to transitions from occupied to unoccupied states. In this way the expression becomes 



Xab 



e h(u} + 27) 



E X/ ' 



Wi-Wf ( (i\v a \f)(f\v b \i) . (i\v b \f)(f\v a \i) 



U>if 



oj + 17 + Wif 



U> + - LOif 



(81) 



D. Application to band models 

To apply this result, we need to use the energy levels 
appropriate for the bands under consideration. The dis- 
cussion is restricted to a parabolic two-band model, with 
the bands separated by some gap energy E g . There are 
effective masses m* h and m* for the occupied (lower) and 
unoccupied (higher) bands, respectively. Each band is 
affected by the Zeeman shift in the same direction; there 
are not Landau level shifts. One can measure energies 
from the top of the lower band. Then the energies Ei for 



the occupied band (E^, valence band) and the energies 
Ef for the unoccupied band (E e , conduction band) are 
assumed to be 

h 2 k 2 

Ei = E h = -—^ - \m % hw B , (82) 
zm h 

Ef = E e = E g + - \m f hLo B , (83) 

(84) 
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These Zccman shifts apply to positive charges: they arc 
reversed in sign for negative charges, taking ujg < 0. 
The azimuthal quantum numbers are m, and m f . They 
are restricted by the orbital angular momentum num- 
bers for each band, U and If, respectively. Assuming 
vertical transitions that conserve linear momentum Sk 
(negligible photon momentum), the transition energies 
are 

ft 2 k 2 

hajif = -Eg - — - - i(mj - m f )hhj B , (85) 
where the reduced mass m* is defined by 

-T = ~T + A ■ (86) 
m m* m h 

We write the transition frequencies in the following man- 
ner: 

oJif = —oj g — s 2 + |Amw B , (87) 

where the gap frequency cj g , scaled wave vector s, and 
change in azimuthal quantum number Am are 



w g = -^-, s=\j^— z k, Am = m f -m t . (88) 

Only momentum-conserving transitions between two 
selected bands at some wave vector k are considered. 
The matrix elements needed are approximated in a form 

hk 

(k'l'm\v x \klm) = — -M(k)<5k',k6i',i±i<W,ra±i • (89) 
m 

The last Kronecker deltas reflect the electric dipole se- 
lection rules, Al = ±1, Am — ±1. The dimensionless 
matrix element M(k) is assumed to be some constant 
for the transitions of interest. 



These velocity matrix elements are proportional to 
corresponding position matrix elements, see (f76j), or even 
the matrix elements of the 7? operator. We only need the 
components of operators along x and y. But the angu- 
lar part of these matrix elements is due to the electric 
dipole selection rules. That angular part has the follow- 
ing symmetries, from matrix elements between spherical 
harmonics, 

(l'm'\vy\lm) = -iAm (l'm'\v x \lm), Am = ±1 . (90) 

This directly affects the susceptibility for each circular 
polarization. From (^0|) we have the diagonal part as 

X*,~5>/*l(/M*>| 2 . (91) 
fi 

but the off-diagonal part as 

Xxy~J2(-iAm)g fi \(f\v x \i)\ 2 . (92) 
fi 

It is clear that x an d e have the same symmetry. Then 
the susceptibilities for the right and left circular polar- 
izations vary like 

XR = Xxx - iXxy ~ E (1 - Am)j/i \(f\v x \i)\ 2 (93) 

fi 

XL = Xxx+iXx y ~^2(l + Am)g fl \(f\v x \i)\ 2 {94) 
ft 

In these expressions, only Am = — 1 (Am = +1) con- 
tributes to xr (xl)- Each factor is a Kronecker delta, 
i.e., (1 ± Am) = 2S m .— mi ±i. The following expressions 
result for integration in the band model expressed us- 
ing the transitions between occupied (lower band) and 
unoccupied (higher band) states: 



2ne 2 



XR 



XL 



2ne 2 



EE 

* / 

o u 

j^a EE 



W i ~ W f |-\|2 

Wif 



f—rrii — 1 



f—mi + l 



ui + i-i + Uif u + ij — Wif 



* f 



U>if 



\(f\ix\i)\' 



Jm f—rrii-\-l 



J-m f —rrii — 1 



(95) 
(96) 



The only difference between these is the swapping of the Kronecker deltas with the helicity index: 
Kronecker deltas. Then the two cases can be written 

o u 



in terms of a single expression, replacing the ±1 in the 2ne 2 . Wi — Wf 2 



w + 27 + usij w + vy — Wif 



(97) 
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To proceed further, it is necessary to evaluate the sums. 
This can be facilitated by converting them to integrals 
over the allowed transitions, which depends slightly on 
the dimensionality of the bands under consideration. 



E. Interband transitions between 
three-dimensional bands 

The band structure of interest could be effectively 
isotropic and three-dimensional, say, for the case of some 
semiconductors near the T point (k = 0). Therefore it 
is interesting to consider the IBT contribution for this 
model, before doing a similar analysis of the reduced 
one-dimensional band model for metals. 

Converting from a sum to an integral with — > 
(2x) 3 / using the assumed form for the matrix ele- 
ments, and then changing to s = \Jhj2m* k as the vari- 
able of integration, the interband susceptibility can be 
written 



Xv 



(98) 



where Q contains all the constant normalization factors, 
and the interband transition integral T v {ui) contains all 
of the frequency and temperature dependence: 



2ne 2 \M\ 2 H V 4tt (2m* 



m 2 o e (2tt) 3 3 \ h 



5/2 



T = 



—EE 



UJif 



Jm e—'rrii—v 



U) + ij + coif uj + 27 — Uif 



(99) 



(100) 



To arrive at this, J dflk 2 = ^-k 2 was used for the an- 
gular part of the int egratio n. The upper limit is a Fermi 
wave vector sf = \J h/2m* kF needed to sum over all the 
states of the occupied initial band. We can let Wi = 1 for 
the lower band, but keep the the temperature-dependent 
occupation probability Wf > for the upper band. In 
this way, any thermal effects due to an initial population 
in the final band will be included. 

The first integral in (100) uses A to = +v and the 
second uses Am = —v. These choices enter in the ex- 
pression for uji / (Am) . In terms of the scaled wave vector 



s or a related excitation variable x — > 



s , one has 



u if (+u) = -u g - s 2 + v— = -x + &, , (101) 

Uif(-v) = _ s2 _ V -7T = ~ x - &> ■ ( 102 ) 



The variable x is the excitation energy above the lower 
band, and £„ = \vuib is the polarization-dependent Zee- 
man splitting. Then the denominators in ( |10C| ) in the 
two terms are (w + i'y + Q v — x) and (w + + C, v + x), 



respectively. This suggests introducing notation for a 
Zeeman-shifted complex optical frequency, for the two 
circular polarizations: 



UJ,, 



: ui + i'y + ( u = ui + ij + 



(103) 



Thus, most of the polarization-dependent effects will be 
carried by this shifted frequency. 

To include th e Fe rmi occupation factor Wi — Wf, the 
sum over mi in ( |100| ) can be done first, holding mj fixed. 
Then both terms will have the same occupation factor, 
taken to be 1—Wf. The first term in ( |100| ) uses only mi = 
mf — v, the second uses only m, = TO/ + v, assuming 
those states exist in the lower band. The final state 
energy, however, depends on the effective mass m* in the 
upper band, whereas the transition energy depends on 
the reduced mass to*. If these are nearly the same, i.e., 
m* w m* h to*, then the final state energy measured 
relative to the top of the lower band can be taken as 

Ej ~ h (u) g + s 2 — |to/w_b) =h(x— gTO/Ws) . (104) 

This will lead to an occupation factor un — w/ for each 
term determined from the Fermi energy Ep, 

g mf {x) w 1 - F(E f - E F ) = 1 - F(x;m f ). (105) 

Inouye et al.0 and also Scaffardi and TochcE! have 
discussed bound electron response in a ID band model 
without B, applying an expression which is an integral 
over the excitation variable x. We can transform the 
integral I u into an integration over x to compare with 
those results. In the simplest case, when the lower band 
has a higher angular momentum Zj than the upper band 
value, If, then both states mi = mf ± v ar e available for 
all the allowed to/. In this case, Eq. ( |100| ) becomes 



T — 



J + 17 



E 

mj 



dx J K i „, , - 9 ' - . (106) 



{x 2 -W B ){x 2 



The function g m . — u>i — Wf includes the magnetic field 
effects on the final state occupation probabilities, which 
are slightly different for each m / level. The upper limit is 
determined by the scaled wave vector at the Fermi level, 
xf = uj g + s 2 F . Although this is the a 3D band model, the 
similarity to the corrected expression from Scaffardi and 
TochotJ is clear. When reduced to one dimension (use 



y/x — uj g in place of (x - 



\3/2 



, together with a different 



constant factor Q out front), it recovers the expression 
from Ref. |l5| for oj b -»■ 0. 

The main effect of the DC magnetic field is to shift the 
optical frequency oppositely for the two circular polar- 
izations. The second important effect is the modification 
of the to/ occupations with B. There appears another 



effect due to the factor x 2 



1 , , 2 

4 U B 



in the denominator, 



however, that is relatively small, quadratic in the field, 
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and not dependent on the polarization. The dependence 
on 2v in the numerator is interesting. 

In the limit where the initial states are fully occupied 
and the final states are fully unoccupied (all Wf = 0), 
these integrals can be carried out analytically, see Ap- 
pendix [A|. In the further approximation where the upper 
limit of integration is set to s F — > oo, the result for this 
complex integral is: 



K-c„) 3/2 -K + C„) 3/2 } 



(107) 



The lower and upper band angular momenta combine to 
give the multiplicity of transitions, g m = li + If. This 
result corresponds to the situation of intrinsic semicon- 
ductor particles where the Fermi level is near or below 
the middle of the gap. This limit takes out any tempera- 
ture dependence. However, that should be a small effect 
only from the rounding of the occupation of levels near 
the Fermi energy. Assuming that the gap frequency uj g 
is large compared to the Zeeman splitting lub , tempera- 
ture effects on the FR could be small. Furthermore, the 
limit of zero damping is simple to read out from this ex- 
pression. However, generally, we can calculate all results 
more precisely using the full theory result in equation 



F. A ID band model 

The band structure in real materials can be very com- 
plicated. The isotropic 3D model just discussed is an 
idealisation for real solids. In the band structure of 
goldES in the L direction at the Fermi surface, the im- 
portant IBTs occur from the d valence band to the sp- 
conduction band, where the two bands are separated by 
a gap of Eg rj 2 eV. Because this happens along a parties 
ular direction, a ID band model is useful. Inouye et a/o 
used a ID band model to get the absorption contribution 
in gold due to the interband transitions for bound elec- 
trons. They were able to fit the absorption very well for 
photon energies above 2 eV up into the ultraviolet, well 
above the small particle plasmon resonance frequency. 
Without interband transitions taken into account, it is 
impossible to get such an accurate description. 

The main difference from the 3D problem, is that the 
integration may start out in 3D, but needs to be reduced 
to an effective integration along only the active direction 
of the band, which is taken as the k x -axis. The details 
of how this is done are not so important. One way is to 
convert the sum in (|97| ) into integration in a box along 
Cartesian axes, using that to define an effective Fermi 
wave vector (slightly different from the standard defini- 



tion) by 



N 



2V 
(2^7 



k f 



2Vk 3 F 



Then an effective Fermi wave vector is defined here as 
k F = 7r(n/2) 1//3 , which differs from the usual expression 
k F = (S^n) 1 / 3 by the factor (tt/6) 1 / 3 « 0.8, of no real 
importance. Now in the sum over wave vectors in (|97j), 
the transverse (inactive) coordinates y, z are integrated 
out, leaving only the integration over the band coordi- 
nate: 



E 



2V 
J2nf 



{2k F f 



dk x — 



Vk 



2 r k 



dk x . (109) 



Normalized this way, this final integration will correctly 
give a dimensionless result for Xu- The conversion from 
sum to integration is carried out by including the factor 
Vk F /ir 3 rather than the more usual factor V/(2ir) 3 used 
for 3D. Then for ID, the interband susceptibility can also 
be written \v = QT v {ui), the latter being an integration 
over s = Jft/2m* k x , and containing the frequency and 
temperature dependence, 
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(111) 



This is very similar to the expression (10C) for 3D. The 
only difference is the factor s 2 for ID in place of s 4 in 
3D. 

The expression for T v can be written with the excita- 
tion variable x = cu g + s 2 ; In the most general case, both 
states nii = rrif ± v do not exist for all choices of m/, 
and the two terms in the integrand cannot be combined 
in a simple form. One must then evaluate (111) by 



T v = 



-1 



1 



^2 dx 9m f (x) 

J UJ n 
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Jrrii—m f —v 



(112) 



When the angular momentum If of the upper band is 
less than that of the lower band, k, both the choices 
rrii = rrif ± v always exist for any m / , and this can be 
written instead as 



dx 



g mf (X) Xy/X - LOg 

(.* 2 - i W |)(z 2 -^ 2 ) 



(113) 



As mentioned earlier, the only change from the 3D ex- 
pression is that the power is now | instead of | in the 
numerator. That is a complex integrand. If the real and 
imaginary parts are needed, it can be written in terms 
of real quantities, better for comparison with the B = 
equations: 



T v = 



u 2v 
to + ij 



^ / dx g mj (x) 



(x 2 + 7 2 - cj|) +2ijuj F 
(x 2 + 7 2 — ujp) 2 + A"f 2 u 2 F 



Xy/x — UJ„ 



(114) 



Here a notation for a real shifted Faraday frequency is 
used, 



UJp = U> + 2 VUJ B , 



(115) 



that is dependent on the polarization index v. In the ap- 
plication of this to gold nanoparticles, we assume transi- 
tions from the lower d-band to the upper p-band, hence 
there are values to/ = — 1,0,1 to be summed over. 

The last result also has an approximate expression 
that applies if the Fermi energy falls well within the band 
gap. Then, the Fermi occupation factors can be approx- 
imated with Wi — Wf — 1, removing all the temperature 
dependence. The integral for T v can be done exactly in 
this case, see the Appendix. There results 



n(k + lf) , , . — 



(116) 



Obviously the expression is similar to that for the 3D 
band model, but now the powers are | instead of |. A 
comparison of the full integral ( |114| ) and this approxi- 
mation (116), as functions of the photon energy, can be 
seen in the Appendix in Figure ^. 
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free electrons: 

ffl p = 1.305 x 10 16 rad/s (X = 144 nm) 
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bound electrons, IBT integral: 
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FIG. 4: (Color online) Fitting of the absorption of 17 nm 
diameter gold particles in water solution, according to the 
ID band model for the interband dielectric response. Fit was 
made via a Monte Carlo search, allowing both the bound 
electron and free electron parameters to be varied. Their 
final adjusted values, for T — 300 K, are indicated on the 
Figure. The fitted gold volume fraction is f s — 5.95 x 10~ r . 
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G. Interband parameters for gold 

Following Inouye et aZ.Q and also Scaffardi and 
TochoEj, the ID band model was applied here for the di- 
electric response of gold nanoparticles. As for the Drude 
model, we can fit the ID interband permittivity, Eq. 
( |114[ ) for B = to the absorption of a solution of 17 
nm diameter gold particles in water. The absorption 
data from 350 nm - 900 nm were fit to Eq. (114) (using 
also the MG theory presented earlier) while allowing the 
gap energy hui g , Fermi energy Ep, bound electron damp- 
ing 7, and normalization constant Q to be varied. We 
also allow the gold volume fraction f s and free electron 
plasma frequency lo p and damping j p to be varied. 

This is a multi-parameter search, which was carried 
out via a Metropolis Monte Carlo algorithm. For the 



FIG. 5: (Color online) The IBT cont ribu tion to the permit- 
tivity, from the ID band model, Eq. (114), using the param- 
eters of Figure ^. In this model both the real and imaginary 
parts remain positive for all frequencies. 



effective energy function to be minimized, it was found 
practical to use the sum of the absolute differences be- 
tween the experimental data at each frequency and 
the theoretical expression a t h, he., J2i \ a i{^i)~ <x*h.( u i)\i 
instead of the squares (this gives a more uniform weight- 
ing but looser fit to the points). A reasonable fit and 
the associated parameters are shown in Figure ^. The 
gap energy E g = 2.02 eV, Fermi energy Ep — 2.53 eV 
and plasma frequency (haj p — 8.59 eV) determined in 
this fit are consistent with the values used in Ref. |l5| . 
The dampings, fvy p — 0.152 eV (7^ = 4.33 fs) and 



ftq = 0.269 eV (7 _1 = 2.45 fs), are somewhat different 
from those for bulk gold, and the resulting real and imag- 
inary parts of dielectric function e(uj) are different from 
those for bulk gold.Eil However, 7 P is consistent with the 
prediction from Eq. ([43]), as expected due to extra sur- 
face scattering and other factors for the nanoscale par- 
ticles. If this surface scattering effect is taken out, the 
model reproduces the real and imaginary parts of e(w) 
for bulk gold that are completely consistent with those 
found by Johnson and ChristyEj 

A value of gold volume fraction f s = 5.95 x 10 -7 was 
needed in the fit. It is about half the value estimated 
by the techniques in Ref. ||(| see Sec. y, showing the 
difficulty in estimation of f s in the lab. We might note 
that this fit is not strongly constrained; it was deter- 
mined only by the absorption data; other values cannot 
be strongly ruled out. Unlike the Drude model presented 
earlier, this model is very good at fitting the ultraviolet 
end of the dielectric properties, and also fits the infrared 
end better. Thus, we expect it should give more reliable 
predictions for the Faraday rotation properties. 



VI. FARADAY ROTATION DUE TO GOLD 
NANOPARTICLES 

Here we apply the ID band model to the interband 
transitions, and compare the theory to experiments for 
gold NPs. Based on the fitting of parameters for 17 
nm diameter gold particles in the previous section, the 
results for scaled complex Faraday rotation angle = 
T + iZ = {ip + iX)/(Bzf s ) can be estimated using Eq. 
( p^[ ) together with the MG theory for the response of a 
dilute solution, Eq. (^). The real and imaginary parts 
are the Verdet and ellipticity factors, respectively, scaled 
by volume fraction f s of gold in the solution. 

The theory results for 17 nm diameter gold NPs are 
compared with experimental data for the FR spectrum 
in Fig. H The experimental data for T — v/f s has been 
scaled down by a factor of 1/10 so that its negative FR 
peak is similar in magnitude to the theory. Although 
some details are not identical between theory and ex- 
periment, the general trend in T as a function of wave- 
length is similar. Both show a strong negative peak due 
to the plasmon near 520 - 530 nm, and a zero cross- 
ing in T near 550 nm. This peak is about ten times 
stronger when using the IBT theory than the simpler 
Drude approach presented in Fig. y, but the negative 
FR peak in the experimental data is still much stronger 
than in both theories. The experimental result for T, 
however, is obtained by dividing raw FR data by the es- 
timated gold volume fraction (using techniques of Ref. 
p6| ), f s = 1.23 x 10~ 6 . Any error in the volume fraction 
will modify the experimental value of T. The volume 
fraction divides out in the theory for T. .iZarious mul- 
tiple scattering and backscattering effects^ and other 
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FIG. 6: (Color online) (a) The real and imaginary parts of 
the scaled complex Faraday rotation $ = (ip + iX)/(Bzf s ), 
from experiment (points), and calculated including IBTs for 
bound electrons in the ID band model for gold particles of 17 
nm diameter (lines) . The real part (solid curve) is the Verdet 
factor per volume fraction, T. The experimental data for T 
has been scaled by a factor of 1/10, which brings its peak to 
the same size as in the IBT theory. The imaginary part is the 
ellipticity factor scaled by volume fraction (not measured). 

similar aggregation effects^ not included in this theory 
could explain this large discrepancy. 

Above 560 nm both theory and experiment indicate a 
positive Faraday rotation angle. According to the theory, 
there is also a strong positive peak in ellipticity expected 
around a wavelength slightly larger than that due to the 
plasmon (around 550 nm). Notably, both the rotation 
and ellipticity tend towards zero at short wavelengths, 
removing the artifact present in the classical Drude ap- 
proach. 



VII. DISCUSSION AND CONCLUSIONS 

The interband electronic transitions are known to have 
a considerable effect on dielectric properties of gold and 
other metals. Although it is popular to consider only a 
simple Drude model for quasi-free electrons, it is shown 
here to be inadequate for describing, say, the absorption 
in gold nanoparticles. A fit for the dielectric parame- 
ters based on a Drude model, combining free electron 
and bound electron contributions, was given in Fig. 
for a dilute solution of gold NPs. The absorption peak 
near 520 nm, associated with the excitation of a plas- 
mon resonance in the NPs, can be fit rather well. In the 
violet and ultraviolet, however, it is impossible to get a 
good match between the theory and experimental data 
for gold particles of average diameter of 17 nm. This ap- 
proach leads (incorrectly) to negative values of the real 
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part of permittivity at shorter wavelengths (below the 
plasmon resonance), which is why the fit fails in that 
wavelength region. 

For the optical effects such as absorption or Faraday 
rotation, it is clear that a correct description of e(w) 
for a macroscopic sample completely determines the out- 
come. If the absorption cannot be fitted properly, then 
it is hard to see how the Drude approach could explain 
Faraday rotation very well. Instead, we include the in- 
terband transitions for the bound electrons, based on the 
quantum formulation due to Boswarva et a/.cJ and also 
Adler.EJ. However, the calculation done here is specifi- 
cally for nanoparticles in the presence of a DC magnetic 
field B. Bulk calculationsEH have assumed Landau lev- 
els due to B for an infinite sample, but the Landau radius 

r ° = f° r wea, k B will be larger than the typical NP 

size. Hence, the Landau levels have no physical sense for 
weak magnetic field; they are not the quantum energy 
eigenstates. The states are modified due to the geometric 
confinement of a limited size nanoparticle. The shift in 
energy states used here is a Zeeman shift due to the DC 
magnetic field for electrons in bands of specified angular 
momentum. 

We are applying this calculation at a magnetic field 
strength such that the Landau radius is only about twice 
the NP radius. This more intermediate field is needed 
to insure adequate signal to noise ratio in the FR data. 
This should be within the range of the quantum theory. 
However, at a magnetic field strength about four times 
larger, the Landau radius will match the NP radius, and 
crossover to a more bulk-like behavior could be expected. 

We have sketched out the details of the IBT integral, 
especially for the ID band model. While it was required 
to fit the amplitude factor Q and other dielectric param- 
eters, this results in a much better fit between theory 
and experiment for the absorption, Fig. [I| The fitting 
leads to an interband permittivity whose real and imagi- 
nary parts remain positive for A > 200 nm, Fig. | Thus 
we expect this approach to be good when applied to the 
Faraday rotation properties of a dilute solution of NPs. 

When applied to a solution of 17 nm diameter gold 
NPs, one finds a strong negative peak in the Verdet fac- 
tor, as expected, near the plasmon wavelength of about 
525 nm. This negative peak is certainly due to extra 
Faraday rotation that is associated with the surface plas- 
mon mode in nanoparticles. This is at least partially 
confirmed by experiments at a small volume fraction of 
gold. The theory also predicts a positive peak in ellip- 
ticity at a wavelength somewhat larger than the plas- 
mon wavelength (nearly 550 nm for 17 nm gold NPs). 
However, at present, the experimental results for scaled 
Verdet factor T = v/f s seem to be about an order of 
magnitude stronger than that indicated by this theory. 
Some of this discrepancy could be due to uncertainty 
in the true value of gold volume fraction. There may 
also be other processes present in the NP data, such as 



aggregation^ and backscattering effects, til that are not 
included in the model. At present we do not have com- 
parable experimental results on the ellipticity factor. It 
will be important to clarify this discrepancy in the size 
of the magneto-optical responses. 
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Appendix A: Evaluation of integrals over scaled 
wave vector 

This analysis applies when the Fermi level is well 
within the band gap, and then the upper band is as- 
sumed to be unoccupied. This removes the temperature 
dependence from the model. 



1. 3D band model 

The integral for the 3D band model, equation (|l00| ), 
is best evaluated using s = ^Jhjlm* k as the variable of 
integration. We assume Wi = 1, wj = 0, and initially 
include the damping parameter 7. For the sum over m/, 
the Kronecker deltas select m f — rrii ± v for each inte- 
grand. For the sum over m„ between bands with defined 
orbital angular momentum li and If, there are multiple 
equivalent transitions, all with the same thermodynamic 
weight. This leads to a multiplicity 

g m = 1 = min(2/ 4 + 1, 2l f + 1) = h + If . (Al) 

mi 

Essentially, this constant factor replaces the temperature 
dependence that would have been included by the occu- 
pation function, g m , (x) in the full theory expression. 
Then in this approximation, the interband susceptibility 
is Xv = QT„(u!), with Q defined earlier in (|99|), and the 
transition integral expressed 

T v = {K x +K 2 ). (A2) 

u + 27 
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T v has been split into two similar integrals, 

[ SF +ds s 4 

'o 



K, = 



(S 2 + UJg ~ C)(S 2 + UJg 

-dss 4 



, (A3) 
(A4) 



l (S 2 + UJg + (v)(s 2 + UJg + W„) 

The shifted frequency uj u is defined in (|l03| ) and C, u — 
\vujb- Consider K\ using the partial fraction expansion, 



ds s 

u v — Cv 



1 



1 



+ UJg — U) V 



S 2 + UJg - C,V 



(A5) 

The integral K 2 is the negative of this with £„ and w„ 
reversed in sign: 



Kn = 



S F dss 4 

u v — C,v 



1 



1 



S 2 +UJg+( l/ 



(A6) 

Note that ui v — (, v = u> + 27. These integrals can all be 
found from the indefinite integral, 



ds s 



1 



+ a a tan" 1 (-) , (A7) 



where the parameter a is complex for the four differ- 
ent cases where this is used. Thus the function /(s) 
is defined with the analytic continuation of the inverse 
tangent to complex arguments. Applying this gives 
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-,„.W 2 



(UJg - UJ V ) 



Kn = 



uj + 17 
(ujg -C,) 3/2 tan 
1 




(ujg + cj„) 3/2 tan 



uj + 17 

(ujg + Cf^t&n- 1 



It is interesting that once these are summed to produce 
T v = J™. {K\ + K 2 ), all linear and cubic terms in s 
cancel out, leaving only the inverse tangents. These are 
evaluated at the upper limit sp. Due to these cancel- 
lations, the upper limit can be let to go to infinity as a 
reasonable approximation (also, the choice of a finite sf 
may be difficult in any case). With s real, the limit of 
the complex inverse tangent can be shown to have the 
following dependence on the parameter a: 

lim tan -1 - = sgn(Re{a})^, (A10) 

lim tan -1 — = sgn(Im{a})^. (All) 

The latter form is useful for the first tan -1 in K\, chang- 
ing y/ujg — uj v to iy/uju — UJg, which is more convenient 
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FIG. 7: (Color online) The IBT contribution to the permit- 
tivity for the ID band model, as functions of photon energy 
in eV, using the parameters of Figure ^[ The solid curves 
apply the full theory, Eq. (114), including the Fermi occupa- 



tion factors, for temperature 300 K. Note the peak in Re{e} 
at the Fermi energy (2.52 eV) in the full theory. The dot- 
ted curves show the result of assuming a total ly un occupied 
upper band, Wf = 0, obtaining e(uj) from Eq. (A21). 



when assuming above the gap excitation. Then the re- 
sulting interband integral for the 3D band model under 
these approximations is 



{i(U V - UJgf /2 + {Uj y + UJg) 3/2 



2(cj + z 7 ) 2 

-K- C ,) 3/2 -K + C,) 3/2 } 



(A12) 



These complex square roots are defined by the root with 
the positive real part, i.e., the root in the first or fourth 
quadrant. It is interesting to notice that the photon fre- 
quency does not appear inside the last factors involving 
the Zeeman shift £„. Also, in the limit of zero damping 
7 — > 0, the very first term is the entire imaginary part. 
For below the gap excitation (w„ < uj g ) the imaginary 
part becomes zero in the absence of damping. For typi- 
cal material parameters, this approximation gives r esult s 
very close to the full theoretical result from Eq. ( |106| ). 
This is mostly due to the high dimensionality; in the ID 
band model this approximation is farther from the full 
theory expression. 



2. ID band model 

The analysis is nearly the same, with interband sus- 
ceptibility expressed via Xv = QT u (uj). One has the 
same expression for the multiplicity, g rn = k + If. The 
weighting factor Q is now given in equation (110). The 



18 



expression for the transition integral can still be written 

9m 



T = 



Ul + 



{Ki + K 2 ) , 



(A13) 



except that in this case, the parts have the power s 2 in 
place of s 4 in their numerators: 



K x = 



+ds s 2 



(S 2 + UJg ~ Cv){s 2 + - LO v ) 



, (A14) 
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The partial fraction expansions are 
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K t = 



ds s 
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(A16) 

and for K 2 , reverse the overall sign and the signs on u> v 
and Cw, 



K 2 = 



1 



1 



/0 u f ~ Cv 

The basic integral needed is 
ds s 2 



S 2 + UJg + U) V S 2 + UJg + C, V 



s — a tan 



(A17) 



(A18) 



Applying this to all the sub- integrals, the results are sim- 
ilar to those in ID, 



Ki 
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V UJ 9+' 



+ (cj s + ^) 1/2 tan- 1 



+ C/ 



(A20) 



Again, it will be useful to reverse the order in the radical 
in the first term in K x , assuming above the gap excita- 
tion. Then we apply again, yjuj g — ui u — > iy/ui v — ui g . 
Letting the upper limit of integration sf — > 00, and in- 
serting the limiting values of the inverse tangents, there 
results 



2(uj + i 1 ) 2 

+ V UJ 9 - + V UJ 9+ C^} • 



(A21) 



Note that for typical values of the parameters, both the 
real and imaginary parts of T v (uj) are positive. This 
function derived has a peak in its real part, for frequency 
near the gap frequency ui g . Othe rwise , it is very similar 
to the more complete theory of (114) that includes the 
varying Fermi occupation factor. That complete theory 
differs primarily in that the location of the peak in its 
real part is near the Fermi energy rather than the gap 
energy. 
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